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We study s-wave pion-pion (tttt) scattering length in lattice QCD for pion masses ranging from 
330 MeV to 466 MeV. In the "Asqtad" improved staggered fermion formulation, we calculate the 
7T7T four-point functions for isospin 1 = and 2 channels, and use chiral perturbation theory at 
next-to-leading order to extrapolate our simulation results. Extrapolating to the physical pion mass 
gives the scattering lengths as m^a 1 ^ 2 — —0.0416(2) and m^a 1 ^ — 0.186(2) for isospin 1 = 2 and 
channels, respectively. Our lattice simulation for tttt scattering length in the 1 = channel is 
an exploratory study, where we include the disconnected contribution, and our preliminary result 
is near to its experimental value. These simulations are performed with MILC 2 + 1 flavor gauge 
configurations at lattice spacing a ~ 0.15 fm. 
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I. INTRODUCTION 

Pion-pion scattering at low energies is elemental and 
important hadron-hadron scattering process. The s-wave 
7T7r scattering lengths are predicted at leading order (LO) 
in chiral perturbation theory (%PT) by Weinberg [l| in 
terms of pion mass, m,, and pion decay constant, f w , as 

m,a[=f » -^-p = 0.160; mW= 2 ~ - -0.0456. 

The next-to-leading order (NLO) corrections depend on 
unknown low energy constants, which can be obtained 
from experimental measurements or lattice QCD. 

The recent measurements of the K e 4 decays |2| and 
K r± -+ 7r ± 7r°7r° decays by NA48/2 at CERN g|] give, 
m n ai=° = 0.221(5) and m^ai= 2 = -0.0429(47). Includ- 
ing xPT constraints in their analysis, the determination 
of s-wave 7T7t scattering lengths reaches [5j, [6j : 

mW=° = 0.220(5); mW= 2 = -0.0444(10). (1) 

Lattice calculations of tttt scattering have been stud- 
ied in quenched QCD by various groups 043 , and first 
full QCD calculation of tttt scattering length was done to 
study isospin 1 = 2 s-wave scattering [13j. First fully- 
dynamical calculation of / = 2 tttt scattering length was 
performed by NPLQCD [H, Mixed-action %PT at 
NLO was used to perform the chiral and continuum ex- 
trapolations, and obtain 

m n ai= 2 = -0.04330(42) and l{= 2 {p) = 6.2 ± 1.2, 

where ^ 2 (/i) is a low energy constant (LEC), which is 
evaluated at physical pion decay constant f 7r _ p hy Using 
the Nf = 2 maximally twisted mass fermion ensembles, 
Xu Feng et al [l6[ adopt the lightest pion mass, perform 
an explicit check for large lattice artifacts, and find 

mj 1 ^ 2 = -0.04385(28), li= 2 (ti = /, iPhy ) = 4.65(85), 

which is in agreement with the above experimental mea- 
surements and phcnomcnological analysis. 



So far, only few efforts have been taken for / = case. 
Using the quenched approximation Y. Kuramashi et al. 
explored the 1 = channel, but the disconnect diagram 
was neglected for some reasons Q. Qi Liu performed 
full QCD calculation for the 1 = channel with the con- 
sideration of disconnected graph, however the scattering 
length has a large error, and serves only as a bound on 
the magnitude [17l |. 

It is well-known that 1 = channel is a great challeng- 
ing phcnomenologically because of a resonance. Encour- 
aged by our trustful measurement of ttK scattering [l8| , 
here we use MILC gauge configurations generated with 
2+1 flavors of Asqtad improved staggered dynamical sea 
quarks [l9| to investigate s-wave tttt scattering for 1 = 
and 2 channels. We measure all the diagrams with extra 
efforts to the disconnect diagram. We note an attrac- 
tive signal for the 1 = channel and repulsive one in the 
1 = 2 case. As presented later, after chiral extrapolation, 
we find at the physical pion mass 

m^ai=f = 0.186(2); m^ 2 = -0.0416(2), 

which is in reasonable agreement with above experimen- 
tal measurements and phcnomcnological analysis as well 
as previous lattice calculations. Moreover, we perform an 
exploratory work for calculating l^ (fi), which is a LEC 
appearing in xPT description of the quark mass depen- 
dence of the scattering length for the 1 = channel. 



II. METHOD 

Let us study the scattering of two Nambu-Goldstone 
pions in the Asqtad-improved staggered dynamical 
fermion formalism at zero momentum. Here we follow 
the original conventions and notations in Rcfs. fjj M. fio} . 
Using operators O v (xi), 7r {x2) for pions at points X\ 
and X2, respectively, we then express the 7T7t four-point 
functions as 

C 7Tn {x 4 ,x 3 ,X2,xi) = (O n (x i )O n (x 3 )Ol(x 2 )Ol r (xi)) , 
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where (• • ■ ) stands for the expectation value of path in- 
tegral. After summing over the spatial coordinates, we 
achieve 7T7r four-point function in zero-momentum state, 

CW(*4)*3)*2j*i) = ^^^^C 7r7r (a)4,a;3,a;2,xi), 

XI X 2 X 3 X 4 

where x x = (xi,ti), x 2 = (x 2 ,t 2 ), £3 = (x3,*3), 2:4 = 
(x4,f4), and t represents time difference, i.e., t = t 3 — t\. 
We choose t\ = 0, t<z = l,t 3 = t, and £4 = t + 1 to avert 
the Fierz rearrangement of the quark lines, and build irir 
operators for isospin 1 = and 2 eigenstates as 

Oi=°(t) = -L{n-(t)Tr+(t + l) + n+(t)n-(t+l)- 

7r°(t)T°(*+l)}. 
0£= 2 (i) = n + (t)n + (t + l), (2) 

with the pion interpolating field operators 

7T+(t) = -^d(x,i) 75 u(x,i), 

X 

7r _ (i) = 2ju(x,t)7 5 d(x,i), 

X 

7r°(t) = =^[u(x,t)7 5 u(x,t) - u(x,t)j 5 d(x,t)}. 

* ^ x 

In the isospin limit, four diagrams contribute to 7T7t 
scattering amplitudes @, H, [HI ■ The quark line diagrams 
contributing to tttt four-point function are displayed in 



Cir7r(*4) *2, t±) 
C^7r(*4, ^3j *2,il) 
C^tt(*4j *3j ^2,il) 



where daggers is the conjugation by the even-odd parity 
(— 1)" for the Kogut-Susskind quark action. We utilize 
the hermiticity properties of the quark propagator G to 
remove the factors of 7 s . The vacuum diagram here is 
accompanied by a vacuum subtraction (2p| . As it is dis- 
cussed in Refs. @, the rectangular and vacuum dia- 
grams create gauge- variant noise, which are suppressed 
by performing the gauge field average without gauge fix- 
ing || [l(| in the current study. 

All the four diagrams in Figure Q] can be combined to 



Figure [TJ and labeled as direct(£>), crossed (C), rectan- 
gular (i?), and vacuum (V) diagrams, respectively. 




(a) D (b) C (c) R (d) V 



FIG. 1: Diagrams contributing to irn four-point functions. 
Short bars stand for wall sources. Open circles are sinks for 
local pion operator. 

It is very difficult to reliably measure the rectangular 
(R) and vacuum diagrams (V) In this work we use 

moving wall source method to solve it @, [To|, HH , namely, 
each propagator, which corresponds to a wall source at 
time slice t = 0, • ■ ■ , T — 1, arc defined by 

J2 D n',n"G t (n")=J2^,{x,t), 0<t<T-l, (3) 

n" x 

where D is the quark matrix. The combination of Gt(n) 
that we apply for 7T7t four-point functions is schematically 
shown in Figure[TJ and we can also express them in terms 
of the quark propagators G, namely, 



(4) 



build physical correlation functions for 7T7t states with 
definite isospin. If we consider that u and d quarks have 
the same mass, the 7T7t correlation function for 1 = and 
1 = 2 can be written in terms of these diagrams, namely, 

Ci=\t) = <0£°(t)|0£°(O)) 

= D+^-C-3N f R+^V, 

Cif(t) = (Oi= 2 (t)\Oi= 2 (0)) = D — N f C, (5) 

where Nf is inserted to address the flavor degrees of free- 
dom of the Kogut-Susskind staggered fermion pj. 



I 

= EE< Rc Tr t G l *3)G tl (a*, h)G\ 2 (an, t 4 )G ta (a*, U)]), 

X 3 X4 

= J2Y.( RC TT ^ G l (^3, *3)G t2 (X 3 , i 3 )<4 (Z 4 , U)G tl (a*, U)]), 

X 3 X4 

= E E< Rc ^[^1 lP*M)G u (aa, h)Gl (as, t a )G tl (x 3 ,t 3 )}), 

X2 X 3 
X2 X3 

Re(Tr[Gf 1 (x 2 , t 2 )G tl (a*, t 2 )])(Tr[G t t 4 (x 3 , t 3 )G u (a&, t 3 )]>}, 
I 
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The s-wave mr scattering length in the continuum is 
defined by 



tan<5o(fc) 
ao = hm ; . 

fc^o k 



(6) 



k is the magnitude of the center-of-mass scattering 
momentum related to the total energy by E^ = 
2 \J m 2 + k 2 of the 7T7r system in a cubic box of size L with 
isospin I. 5o(k) is s-wave scattering phase shift, which 
can be calculated by the Liischcr's formula [22I l23j| . 



/ tan<5o(fc) 



4n f k 2 

■ ^oo J-, ■ 



(7) 



irL -u " V" (2tt/L) 2 
where the zeta function Zqq(1; q 2 ) is defined by 



here q = kL /(2ir), and Z 00 (l;q 2 ) can be efficiently com- 
puted the method discussed in Refs. (13, Il8j. 

The energy can be obtained from the irir four- 
point function denoted in Eq. (O with the large t. At 
large t this correlator will fall as 



C'L(t) oc e 



-E ww t 



(9) 



where E n7r is the energy of the lightest two pion state. 
In the usual manner, pion mass can be evaluated 
through 



(10) 



In our concrete calculation we evaluate the energy shift 
SE = E — 2to t from the ratio 



XM _ c*(o,i,m + i) 



C^(0,t)C^(l,t + l) 



X=D,C,R, and V, 



where C w (0,t) and C ff (l,t + 1) are the pion two-point 
functions. The amplitudes which project out the 1 = 
and 2 isospin eigenstates can be written as 

Ri=o(t)=R D (t) + ^R c (t) - 3NfR R (t) + ^ y (t), 
Ri=2{t)=R D {t)-N f R c {t). (11) 

We extract the energy shift SE from the ratio 

R I (t) = Z I e- SEt + --- , (12) 

where Zj stands for wave function factor || EH • 

III. SIMULATION RESULTS 

We use the MILC lattices in the presence of the 2 + 1 
dynamical flavors of the Asqtad-improved staggered dy- 
namical fermions, the description of its simulation pa- 
rameters are given in Refs. [19l [24[. We analyzed the 



7T7T four-point functions on the 0.15 fm MILC ensemble 
of 360 16 3 x 48 gauge configurations with bare quark 
masses am u d/am s = 0.0097/0.0484 and bare gauge cou- 
pling 10/g 2 = 6.572, which has an inverse lattice spacing 
a- 1 = 1.358^2 GeV [Hill. The masses of the u and d 
quarks are degenerate. 

We use the standard conjugate gradient method to ob- 
tain the required matrix element of the inverse fcrmion 
matrix. We compute the propagators on all the time 
slices t = 0, • • • , T — 1 of both source and sink. After 
averaging the correlators over all T = 48 possible values, 
the statistics are significantly improved because we can 
place the pion source at all time slices. 

With same configurations we compute the irir four 
point correlation functions using six u valence quarks, 
namely, am x = 0.0097, 0.01067, 0.01261, 0.01358, 
0.01455, and 0.0194, where m x is the light valence u 
quark mass. 

In Figure [2] the individual ratios, which are defined in 
Eq. (fTTjl corresponding to the diagrams in Figure HJ R x 
(X = D,C,R and V) are displayed as functions of t for 
am x = 0.0097. It is extremely noisy for the disconnected 
diagram(y), but still we can get a good signal up to time 
separation t = 14. Clear signals observed up to t = 19 
for the rectangular amplitude and up until t = 14 for 
the vacuum amplitude demonstrate that the method of 
wall source without gauge fixing used here is practically 
applicable. 

The values of the direct amplitude R D is quite close to 
unity, indicating a weak interaction in this channel. The 
crossed amplitude, on the other hand, increases linearly 
up to t ~ 17, implying a repulsion in the 1 = 2 channel. 
After a beginning increase up to t ~ 4, the rectangular 
amplitude shows a linear decrease up until t ~ 20, sug- 
gesting an attractive force between two pions. Moreover, 
the magnitude of this slope is analogous to that of the 
cross amplitude but with different sign. These character- 
istics are what we want 0, 0] ■ 

The vacuum amplitude is negligibly small up to t ~ 
10 — 14, and loss of signals after that. This character- 
istic is in well accordance with the Okubo-Zweig-Iizuka 
(OZI) rule and xPT in leading order, which expects the 
disappearing of the vacuum amplitude 0]. Moreover, its 
the errors should be approximately independent of time 
separation t, and increases exponentially like e 2mwt [||. 
Therefore, it is very difficult to obtain its proper infor- 
mation from large time separation. 

In Figure [3] we plot the ratio Ri(r) projected onto the 
isospin 1 = and 2 channels for am x = 0.0097, which are 
denoted in Eq. (fTTj) . A decrease of the ratio of Ri—2(t) 
indicates a positive energy shift and hence a repulsive 
interaction for the 1 = 2 channel, and an increase of 
Ri=o(t) suggests an attraction for the 1 = case. A dip 
at t = 3 for the 1 = channel can be clearly noted (Io| . 

In this work, we use Eq. (fT2j) to extract the energy 
shift SEj, and then insert them into the Eqs. (0) and (0) 
to obtain the scattering lengths. In this work, the energy 
shifts aSE are picked up from the effective energy shift 
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FIG. 2: Individual amplitude ratios R x (t) as functions of t. 
(a) Direct diagram shifted by 0.8 (diamonds) and vacuum di- 
agram (octagons); (b) crossed (squares) and rectangular (di- 
amonds) diagrams. 



plots, and they were selected by searching for a "plateau" 
in the energy shift as a function of the minimum distance 
D min as well as a good fit quality (namely, x 2 ) [HI- 

We utilize the exponential physical fitting model in 
Eq. (|T2|) to extract the desired energy shifts for both 
1 = 2 and channels. In Figure [3] we display the ra- 
tio Ri(t) projected onto both channels for am x = 0.0097. 
The fitted values of the energy shifts, SEj in lattice units, 
fitting range, and wave function factor Zj arc summa- 
rized in Table |U The third block shows energy shifts in 
lattice units, Column four shows wave function factor Z/, 
Column five shows time range for the chosen fit, and Col- 
umn six shows degrees of freedom (dof) for the fit. The 
wave function Zj factors are pretty close to unity and 
the x 2 /dof is pretty small for the 1 = 2 channel, indi- 
cating the values of the extracted scattering lengths are 
substantially reliable, and the Z factors are also near to 
unity, and the x 2 /dof is reasonable for the 1 = channel, 
suggesting the value of the extracted scattering lengths 
are enough safe. 

In our previous work (Lsl . |25| , we have measured the 
pion masses (m^) and the pion decay constants / w [26j . 
which are summarized in Table [TTJ The second and third 
blocks show pion masses in lattice unit and in GeV, re- 
spectively, and Column four shows the pion decay con- 




FIG. 3: Ri(i) for tttt four-point function calculated without 
gauge fixing for am x = 0.0097. Solid line in I = 2 is exponen- 
tial fits for 10 < t < 15 and Solid line in / = is exponential 
fits for 7 < t < 14. 



TABLE I: Summary of simulation results for energy shifts. 



I 


m x 


SE 


Zj 


Range 


x 2 / dof 




0.00970 


-0.0241(19) 


0.825(7) 


7- 


14 


12.8/6 




0.01067 


-0.0235(32) 


0.835(23) 


7- 


14 


9.46/6 





0.01261 


-0.0222(29) 


0.854(19) 


6- 


12 


6.94/5 




0.01358 


-0.0223(30) 


0.856(20) 


6- 


12 


5.33/5 




0.01455 


-0.0217(32) 


0.863(20) 


6- 


10 


4.66/3 




0.01940 


-0.0199(30) 


0.884(19) 


6- 


12 


16.5/5 




0.00970 


0.00670(15) 


0.977(2) 


10 


- 15 


1.07/4 




0.01067 


0.00662(15) 


0.980(2) 


10 


- 15 


0.17/4 


2 


0.01261 


0.00648(14) 


0.984(2) 


10 


- 15 


0.05/4 




0.01358 


0.00640(14) 


0.986(2) 


10 


- 15 


0.03/4 




0.01455 


0.00625(13) 


0.986(2) 


10 


- 15 


0.02/4 




0.01940 


0.00594(12) 


0.993(2) 


10 


- 15 


0.08/4 



stants in lattice units. Now we can substitute these cn- 



TABLE II: Summary of the pion mass and the pion decay 
constants. 





arrin 


m 7r (GeV) 




0.00970 


0.2458(2) 


0.334(6) 


0.12136(29) 


0.01067 


0.2575(2) 


0.350(6) 


0.12264(34) 


0.01261 


0.2787(2) 


0.379(7) 


0.12425(27) 


0.01358 


0.2890(2) 


0.392(7) 


0.12482(32) 


0.01455 


0.2987(2) 


0.406(7) 


0.12600(26) 


0.01940 


0.3430(2) 


0.466(8) 


0.12979(27) 



ergy shifts in Table U into the Eq. ([6]) to achieve the 
scattering lengths. The center-of-mass scattering mo- 
mentum k 2 in GeV calculated by E\ K = 5Ei + 2-m^ = 
2y // m 2 r + k 2 and then its s-wave scattering lengths ao ob- 
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TABLE III: Summary of lattice simulation scattering lengths. 



Isospin 


m x 


fc 2 [GeV] 


do 


m^ao 




0.00970 


-0.0107(4) 


2.95(19) 


0.724(48) 




0.01067 


-0.0109(7) 


3.03(32) 


0.781(85) 





0.01261 


-0.0112(7) 


3.16(33) 


0.882(92) 




0.01358 


-0.0117(8) 


3.40(38) 


0.983(107) 




0.01455 


-0.0118(8) 


3.67(38) 


1.095(112) 




0.01940 


-0.0124(9) 


3.78(50) 


1.297(173) 




0.00970 


0.00306(7) 


-0.495(5) 


-0.121(2) 




0.01067 


0.00316(7) 


-0.509(5) 


-0.131(1) 


2 


0.01261 


0.00335(7) 


-0.537(5) 


-0.150(1) 




0.01358 


0.00343(7) 


-0.548(5) 


-0.159(2) 




0.01455 


0.00345(7) 


-0.552(5) 


-0.165(2) 




0.01940 


0.00378(8) 


-0.598(6) 


-0.205(2) 



tained through Eq. ([6]) for both 1 = 0, and 2 channels arc 
summarized in Tabic Hill The errors come from the statis- 
tic errors of the fitted values of the energy shifts. The 
third block shows center-of-mass scattering momentum 
k 2 in GeV, Column four shows s-wave scattering lengths 
in lattice units, and Column five shows pion mass times 
s-wave scattering lengths. 



We adopt the formula predicted by xPT at NLO to 
extrapolate 7T7t scatteri ng l eng ths to the physical point. 
As suggested in Refs. p^ - llq . we carry out the chiral 
extrapolation of m^a^ 2 and m^a^: in terms of m^/ fit 
instead of m T . Thus we use the continuum xPT forms 
of ag =2 and aQ =0 , which are directly constructed from 
Appendix C in Ref. [27|, as 



j=i 



7< 

167T/2 

8^/ 2 



16tT 2 / 2 

ml 
16n 2 f 2 



91n 7T -5-^(A* = A, P hy) 

9 

TfL 

3 In -jf - 1 - ll^ 2 (p = U,phy) 

J7V 



(13) 
(14) 



where we plugged in the values of the pion mass m^, and 
the pion decay constants f n , which are summarized Ta- 
bic HH and l^ (fi) and l^inl are related to the Gasser- 
Lcutwylcr coefficients Zi as [27[ 



40 r 80 r 

— h -I 1 2 

21 1 21 

8 r 16 F 
3 h + Y h ~ 
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h + 4/ 4 + 9 In ■ 



-"%phy 

n 2 



l 3 -4Z 4 + 31n- 



Trophy 
M 2 



,(15) 



(16) 



here Eq. (Tho)) were extensively used in Refs. |14Hl€ 

The chiral extrapolation of irir scattering lengths, 
m-7T<il =2 and m^a^ are plotted by the dotted lines as a 
function of m 2 in Figure 21 The fit parameters l^°(fj,), 
l^ 2 ([i), and the s-wave scattering lengths m^a n at the 
physical points, where we adopt the latest PDG [28| val- 
ues, are also summarized in Table HVl The chiral scale is 
taken as [i = /^phy From Figure |U we can note that 
our lattice simulation results for the 1 = 2 scattering 
length agrees well with the one-loop formula, while scat- 
tering length for I = have a large error, and is also in 
reasonable agreement with xPT at NLO. 



TABLE IV: The fitted m^ao at the physical point. 



Isospin 


X 2 /dof 




rriir do 


1 = 


0.268/5 


18.674 ± 1.213 


0.186(2) 


1 = 2 


9.864/5 


11.587 ±0.871 


-0.0416(2) 




0.00 



0.05 



0.10 0.15 
m^deV 2 ] 



0.20 



0.25 



FIG. 4: m^-dependence of nn scattering lengths m n ao for 
1 = and 2 channels. The dotted lines give x?T predictions 
at NLO. The cyan diamond points indicate its physical values. 



Although the fitted value of l^.~ 2 (/it) is larger than that 
of other lattice studies EMI, 

our fitted value of m w aiz 
is reasonable consistent with other lattice studies [1J- 
[l6j . Since it is an exploratory study for / = chan- 
nel, there are no lattice comparisons with our fitted val- 
ues of and m^ao- Anyway, our fitted value of 
s— wave 7T7T scattering length for 1 = channel is in rea- 
sonable agreement with the experimental measurement 
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in Eq. ©. 

IV. CONCLUSIONS AND OUTLOOKS 

In this work, we performed a lattice study of s-wave 
tttt scattering for isospin 1 = and 2 channels. We 
evaluated all of the four diagrams, and observed an at- 
tractive signal for the I — channel and an repulsive 
one for the 1 = 2 channel, respectively. Extrapolating 
toward the physical point yields ra^a 1 ^ 2 = —0.0416(2) 
and m-Trd^ = 0.186(2) for I — 2 and channels, respec- 
tively, which are in good agreement with %PT at NLO, 
and m^a 1 ^ 2 = —0.04 16( 2) is reasonably consistent with 
other lattice studies 111! 10 . Moreover, we give an ex- 
ploratory fitted value of s— wave tttt scattering length for 
1 = channel, which is in reasonable agreement with the 
recent experimental measurement Q. 

It is quite stimulating that tttt scattering for the 1 = 



channel can be trustworthily calculated by wall sources 
without gauge fixing. It gives us hope that we can use 
this technique to tackle a resonance, which is still poorly 
understood from lattice QCD. 

In our previous works [25], |29| , [301 ] , we studied and eval- 
uated a mass, and found that the a meson is heavier 
than the tttt threshold for enough small u quark mass. 
These works as well as our lattice investigation for the 
7T7r scattering in the 1 = channel will stimulate people 
to investigate the decay mode a — > tttt. 

Since our study is restricted only at zero momenta, 
it can not provide adequate information on a meson. To 
investigate the a resonance, we should study the tttt scat- 
tering length at the 1 = channel with non-zero momen- 
tum, which is an indicative of a a pole. We are begin- 
ning this kind of work, and the measurement of the tttt 
scattering for the 1 = channel with the momentum 
p = (1, 0, 0) is in progress. 
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